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SUMMARY 


In  this  report  an  inviscid,  compressible,  linearized  disso¬ 
ciating  gas  flow  over  slender  bodies  of  general  cross  section  is 
considered.  Essentially  it  is  an  extension  of  the  conventional 
slender  body  theory  include  the  dissociation  effects.  Both  super¬ 
sonic  and  subsonic  cases  are  examined.  To  our  present  approximation 
the  lateral  forces  and  moments  are  found  the  same  as  those  for  the 
conventional  case,  but  there  exist  additional  non -equilibrium  terms 
in  the  drag  expression.  If  the  body  is  pointed  at  both  ends,  the 
non -equilibrium  term  is  formally  the  same  for  both  the  supersonic 
and  the  subsonic  cases  except  that  the  sign  in  the  respective  cases 
is  just  the  opposite.  Physical  aspects  of  the  nonequilibrium  drag 
problem  are  briefly  considered.  These  considerations  show  that  in 
the  subsonic  case  the  nonequilibrium  term  represents  a  drag  and  in 
the  supersonic  case  it  is  a  thrust.  Finally,  the  supersonic  drag 
of  slender  bodies  of  revolution  at  small  angle  of  attack  has  been 
calculated,  including  dissociation  effects. 
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x,y,z 


Velocity  of  sound 

Contour  of  the  cross  section  in  plane  normal 
to  x-axis  (see  below) 

Drag  coefficient 

Drag 

Inverse  Fourier  transform 

Modified  Bessel  function 

Ratio  of  the  equilibrium  speed  of  sound  to  the 
frozen  speed  of  sound 

Length  of  the  body 

Inverse  Laplace  transform 

Frozen  Mach  number 

Normal  to  the  body  contour  in  y-z  plane 

Pressure 

Velocity 

Radial  coordinate 
Body  radius 
Thickness  of  the  body 
Distance  along  C 
Laplace  transform  variable 
Cross  section  area 

Velocity  component  in  x,  r,  direction 

Rectangular  coordinates,  x-axis  being  in  the  free 
stream  direction 
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Free  scream  velocity 
Complex  potential 
defined  in  Eqs.  (3),  (9) 

Angle  of  attack 

Defined  in  Eq.  (3a) 

Euler's  constant 
Complex  variable 
Angular  coordinate 
Mass  density 

Effective  relaxation  time,  ref.  (1) 

Perturbed  potential 

Fourier  transform  of 

Stream  function 

Fourier  transform  variable 
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Equilibrium  state 
Frozen  state 
Undisturbed  state 
Real  part 


Superscript 
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I.  INTRODUCTION 


The  linearized  inviscid  compressible  dissociating  flows 

of  a  diatomic  gas  are  irrotational  and  satisfy  a  generalized 

wave  equation  of  third  order  in  terms  of  the  velocity  potential 

<P  .  This  equation  was  first  derived  by  Moore*-  and  has  the  same 

form  as  the  equations  for  the  propagation  of  elastic  waves  in  a 

material  subject  to  the  relaxation  of  stresses,  as  obtained  by 
2 

Morrison  , 

13  4  5 

By  means  of  this  linearized* equation,  several  authors  *  *  * 

\ 

have* discussed  some  two-dimensional  flow  problems.  Except  for 

3  14  5 

the  case  of  a  flow  past  a  wavy  wall  ,  it  was  difficult  ’  ’  to 
determine  the  flow  field.  Therefore,  in  most  of  the  recent  results, 
only  the  conditions  on  the  body  surface  were  evaluated.  For  the 
study  of  an  entire  flow  field,  Moore  and  Gibson  further  simpli¬ 
fied  the  generalized  wave  equation  to  a  variant  of  the  telegraph 
equation  for  which  solution  may  be  obtained  in  a  closed  form. 

The  present  report  deals  with  linearized  dissociating  flow 
past  slender  pointed  bodies  of  general  cross-section.  Essentially 
it  is  an  extension  of  the  earlier  works  of  Ward  and  Adams  and 
Sears**'^,  to  include  the  dissociation  effects.  Both  the  supersonic 
and  subsonic  cases  have  been  considered.  The  Laplace  transform 
method  is  used  for  the  supersonic  case,  the  Fourier  transform 
methods  for  the  subsonic  case.  It  will  be  seen  that  the  trans¬ 
formed  equations  (transformed  with  respect  to  x)  are  identical  in 


form  as  those  obtained  by  Ward  and  Adams  and  Sears^’^.  Then 
by  the  usual  slender  body  approximations,  the  problem  can  be 
solved  without  any  additional  simplifications. 

Non -equilibrium  processes  lead  to  irreversible  increase 

in  entropy  along  streamlines  in  the  flow.  However,  this  in- 

8 

crease  of  entropy  is  of  second  order  of  magnitude,  consequently, 
to  our  linear  approximation,  the  flow  is  still  isentropic. 

g 

Kusukawa  and  Li  have  shown  that  the  nonequilibrium  drag  due  to 
the  increase  of  entropy  is  of  the  third  order  in  perturbation 
quantities.  Such  a  nonequilibrium  entropy  drag  will  not  be  con- 
sidered  in  the  present  paper. 

We  remark  also  that  though  we  shall  always  refer  our  dis¬ 
cussion  to  dissociation  in  this  report,  actually  all  the  results 
apply  equally  to  the  case  of  vibrational  relaxation  because  the 
linearized  equation  for  these  two  cases  is  of  the  same  form. 

Development  of  the  present  theory  was  completed  in  late 
spring  of  1961.  Essentials  of  this  work  have  been  recorded  pre¬ 
viously  in  Wang's  thesis*-^  which  was  carried  out  under  Li's 
direction.  Independently,  Clarke^  studied  the  present  problem 
and  obtained  results  similar  to  the  present  theory.  However,  in 
Refs.  (11)  and  (12)  the  authors  held  different  opinions  on  the 
matter  of  the  subsonic  nonequilibrium  drag  on  a  doubly  pointed 
slender  body.  Clarke  concluded  that  in  such  a  case  the  subsonic 


3 


nonequilibrium  drag  is  positive.  Wang  contended,  on  the  other 
hand,  that  the  subsonic  nonequilibrium  drag  expression  can  take 
either  positive  or  negative  sign  depending  on  the  body  shape 
(see  Appendix  of  the  present  paper,  also  see  Ref.  11) .  To  be 
sure,  the  subsonic  nonequilibrium  drag  problem  was  studied  by 

g 

Kusukawa  and  Li  ,  using  a  thermodynamic  approach.  These  recent 
calculations  would  support  Clarke's  conclusion.  Publication  of 
the  present  paper  is  undertaken  at  this  time  since  it  would  com¬ 
plement  the  research  results  of  Refs.  8  and  12.  It  is  hoped 
that  a  quantitative  analysis  of  the  nonequilibrium  drag  effects 
will  be  made  available  shortly. 


II.  SOLUTIONS  OF  THE  EQUATION 


Consider  a  cylindrical  coordinate  system  as  shown  in 
Fig.  I,  the  free  stream  Velocity  U  being  aligned  with  the 
x-axis,  while  the  pointed  nose  being  located  at  the  origin. 
The  body  is  slender,  subject  to  the  same  restrictions  as  in 
the  conventional  slender  body  theory**' Its  cross  section 
may  not  be  circular.  The  medium  is  taken  to  be  a  diatomic 
gas  such  as  pure  oxygen  or  nitrogen  which  forms  a  binary  mix¬ 
ture  when  partially  dissociated.  In  a  linear  theory  the  flow 
is  irrotational,  a  perturbation  potential  <p  may  be  defined: 

?  ‘  V/ 


where  is  the  perturbation  velocity.  For  steady  flows, 
satisfies  the  following  equation^- 


where 


+  ^  +i  +i  ^  1 

^  OL* 

•4  ft  4  1  ^  1 


(1) 


f 


u 


_ 


*  C 
v- 


/C  denotes  the  relaxation  time,  the  frozen  Mach  number, 

and  are  respectively  the  equilibrium  speed  of 

sound  and  the  frozen  speed  of  sound.  It  can  be  shown  that  0^^ 
is  greater  than  <X„  so  that,  generally,  k  is  less  than  unity. 

C  tjm 
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For  the  ideal  cases  of  equilibrium  and  frozen  flows,  there  is 
only  one  speed  of  sound,  and  respectively.  Eq.  (1) 

consists  of  two  parts  (in  two  separate  brackets) ,  one  for 
frozen  flow,  the  other  for  equilibrium  flow.  When  the  relaxa¬ 
tion  V  is  long,  the  frozen  part  dominates;  in  the  opposite  case, 
the  equilibrium  part  will  govern  the  motion.  In  the  following, 
we  shall  refer  to  flows  as  being  "supersonic*1  or  "subsonic"  where 
both  (Aj.  and  are  greater  or  less  than  unity. 


(a)  Supersonic  case 

For  the  supersonic  case,  in  front  of  the  frozen  Mach  line 
from  the  nose  the  disturbances  and  the  perturbation  potential 
vanish.  In  the  disturbed  flow  region,  Eq.  (1)  can  be  written  as 


fc*  L  f  a* 


a- a) 


ft*  »>>3[  jMV-J 

Eq.  (1-a)  can  be  solved  by  the  Laplace  transform  method.  With 
the  notation  S  -  ,  the  operational  form  of  Eq.  (1-a)  is 

^  ~ 


^ +  _L  ^  $>  .«  7 


(2) 


<P 


f  -Sa 

C  clx 


where 
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and  we  have  defined  here 


Since  K  4  1,  hence  ;  may  be  positive  or  negative,  but 

always  real.  It  is  interesting  to  see  that  the  transformed 
equation  (2)  assumes  a  form  identical  to  that  for  the  conventional 
case  of  non-dissociating  flow.  The  difference  appears  only  in  . 
For  purely  frozen  or  equilibrium  flows  Eq.  (3)  reduces  to 

or 

-0 

The  solution  of  Eq.  (2)  may  be  readily  obtained,  for  ex¬ 
ample,  by  the  separation  of  variables  method  in  terms  of  modified 
Bessel  functions  and  (j  •  To  have  an  outgoing 

wave  propagating  downstream  from  the  body,  we  take 

gov*©* 

where  the  coefficients  and  are  to  be  determined  presently. 
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Mow  we  introduce  the  slender  body  approximation.  For 
slender  bodies  -£  ,  where  i  is  the  length  of  the  body. 

Hence  to  describe  the  flow  on  and  near  the  body,  we  retain 
only  the  leading  term  in  the  series  expansion  of  ,  i.e., 


where 


K„(ui>0  4  +*) 


r  =■ 


CKlA)ft 

.5772...  is  the  Euler's  constant. 


n  *  1,  2,... 


Thus 


\oi\X 


<p  =  log  2 

■t  2.  2.^  1  C'n  » 

-r\  -  l 

By  inverse  transform 

<f>  =  <^0  (x  )  lea  -V  C>,  (*)■+• 

O'  T-I 

where 


C^(s)cr>\ir'  0  ^  3>r(s)  SW\  Yg 


Ck)  x ) 


r\ 


(ki  *V'  J 


(4) 


3^(0  Vv,- 

— ■  i — -  -V  _ 


7l 


(5) 


(5a) 


■'a  CO-  „f[-4>Cs>] 

«■<>(*)  =  }''[-ce<s>Os  Hi  +f). 

=  J''[  t*-4  0o->M 

(O  =  j-|  ^ (\,.o ; 

Introducing  <3  =  4  v  z  and  a  complex  potential  W(  ^  ,  x),  we 

have 


=f>JMfr  *©C<>k<^ 

<5 


Ti  i  (6) 


where 
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Q-n CO  —  A ^  L  3  y\  C  *0 


(6a) 


^oCx>;^(x)wiii  be  determined  by  the  boundary  conditions.  In 
particular,  (pc.)  can  be  determined  in  exactly  the  same  way 
as  in  Ref.  (6),  so  that 


0*0  - 


ia^,  s' CO 


air 


(7  a) 


The  term  b„  (x)  can  be  evaluated  in  terms  of  a  (x) .  From 


Eq.  (3),  we  have 


t< 


ki 


5 


,  ion  ujiis 

cT  *  K  &  ±  0  e  +  &. 


2  S  ’3  "  1  o  L£+?» 

The  inversion  of  these  terms  are  found^  to  be 

r'(s)  --  s'co- 

-cH-o 


-  e 


where  C(x)  is  the  Dirac  Delta  function.  Then  by  convolution, 
we  obtain  from  Eq.  (5a) 

G-oC*)  =  a.o0>  tog  Cc-3  • 


+  2  f*oG0 


-?,C*  **0 


(7b) 


.=j& 


X  -  ^ 


The  term  involving  a(0)  has  been  omitted  because  a^CoS^  sfeao 
for  pointed  nosed  bodies.  The  first  two  terms  are  same  as  those 
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given  by  Ward**  except  for  the  appearance  of  the  frozen  Mach 
number  In  the  present  result.  The  last  term  gives  the  non- 
equilibrium  effects.  Any  calculation  not  involving  bD(x)  will 
be  free  from  the  effects  due  to  non -equilibrium  dissociation 
because  it  is  the  only  place  the  non-equilibrium  effects  enter. 


(b)  Subsonic  case 

For  the  subsonic  case,  all  disturbances  die  out  at  a  large 
distance  from  the  body,  certainly  they  vanish  at  x  -  ±  . 

Eq.  (1)  can  be  solved  by  the  Fourier  transform  method.  If  £ 
denotes  the  Fourier  transform  of  <f)  with  respect  to  x,  i.e., 

$  =  f  doc  - 

•  -v 

then  the  transformed  form  of  Eq.  (1)  becomes 


ti  a  i  ±£  4  ±  v2-! 


(8) 


where 


with  p,  and  similarily  defined  as  in  Eq?  (3).  Eq.  (8) 
differs  from  that  considered  by  Adams  and  Sears^  only  in  , 
which  now  is  a  complex  function  of  £0  .  This  equation  can  be 
solved,  under  the  restriction  >  ©  ,  in  exactly  the  same  way 

as  in  the  supersonic  case.  Thus  we  obtain 

$  ~  (*<■*)[  c^C^o)  4.  Sw  to  e  J 


The  modified  Bessel  functions  with  complex  argument  have 
similar  asymptotic  forms  to  those  with  real  argument.  For 
slender  bodies  we  retain  only  the  first  term  in  the  expansion 
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series  of  \<^  ,  the  solution  becomes  formally  identical  with 
that  for  supersonic  flow  given  in  Eq.  (4) . 

(10) 


Making  inverse  transform; 


<KX>  *  =  CO  U  oC*  )  ^  ( 


V 


(11) 


where 


- 

loM  z 


\  Cx)  - 


f*'L-  c°H 

F~' [c„(o)(l°3  |  +  <rj] 


(11a) 


B-„  W  *  F'T  2""‘  Cr 


The  inverse  Fourier  transform  of  a  function  is  defined  by 

r-'M  *  ^  ( -f- o>  e  ’  ,v< 

~  V* 

A^x)  and  BjjCx)  are  real  because  <j>  is  real.  Again  let  ^  -  >,& 
and  consider  <f>(oc,>llb)a8  the  real  part  of  a  complex  potential 
W(  5  ,  a)  so  that 

V  Q  _  /V  'i 

(12) 


Vl(s,.x)--  ^«.0Cx)+2. 

■'  1V»I  - 


1 


Is 


{  l 
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where 


(K-h  00  =  A-n  GO  0=0 

&0GO  ^v\0(  )  mu8t  be  determined  as  before  by  the  boundary 
conditions,  while  Qr0C*S)  may  be  evaluated  in  terms  of  (x)  . 

To  evaluate  ft-0  (x)  ,  we  rewrite  Eq.  (9) 


4  -  cS0-«?) 


*  x  r  co4 

£0*  -V 


'4 


i© 


e  *  0  =  Ccw'1  Cgr^JfQ 


4  Co 


for  >  0  *  we  Put 

•JL 


^  =  I  to  |  /T^ 

so  that 


CO 
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'/4 
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m  -  + *°g  +  ^ 


4  l 


"U  <C  G^I 

9 
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(13) 


The  inverse  transform  of-  these  terms  can  be  shown’*  10  to  be 


F"£fc3H  =  -■ ^  ' 

F"’[(cS  ^ZlJ  =  £CO  ioc 
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rr 


iy 


CO*M  ^  -I 
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-?i 


(14) 


rJ[  ^ 


.  c, 


dT' 


-4  I* 
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~  9silal  —  ^3(|- 
-  C 
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where  £(x)  is  the  Dirac  Delta  function  and  S^t^is  defined 
by 
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sgn  x 


1 

-1 

.  0 


x  >  0 
x  ^  0 
x  =  0 


then  by  convolution,  we  finally  obtain 

Q-o  C*-)  -  ^  -  °»C<) 

2  /? 


) 


I 


2  [  ^  W  ~  J  ^  (3-x^  ^ 


■+1  *oC*) 


e. 


cA<^ 


0<3C  <  i 


(15) 


To  save  writing,  we  have  heare  used  the  condition  that  aD(x), 
a0'(x)  vanish  for  (,<£  x  and  x  <:  0  as  a  direct  result  of 
Eq.  (7a).  Except  for  the  appearance  of  the  frozen  Mach  number 
the  first  five  terms  are  the  same  as  those  given  by  Adams  and 
Sears^.  Again  the  last  term  of  bQ(x)  gives  the  nonequilibrium 
effect.  We  observe  that  the  nonequilibrium  term  has  the  same 
form  for  both  the  supersonic  and  subsonic  flows. 


From  the  classical  hydro dynamical  theory,  we  recognize 
immediately  that  the  first  term  in  Eqs.  (6)  and  (12)  represents 
a  two-dimensional  source  while  the  summation  terms  represent  a 
two-dimensional  doublet  (n  ■  1)  and  higher  order  signularities 
(n  >1).  Therefore  W-  &0C*-)*ePresents  an  incompressible  two- 
dimensional  potential  flow  past  a  cylinder  of  general  cross 
section.  . 
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The  none qui librium  effect  appears  only  in  the  last  term 
of  b0(x).  Since  Eqs.  (6)  and  (12)  differ  from  the  corres¬ 
ponding  formulus  in  the  conventional  slender  body  theory  only 
in  b0(x),  it  follows  that  the  nonequilibrium  dissociation 
affects  only  the  x-component  of  velocity  u',  not  v'  or  w'  in 
the  radial  or  angular  direction.  An  immediate  consequence  of 
this  is  that  the  nonequilibrium  dissociation  does  not  change 
the  lateral  forces  and  moments  because  they  depend  on  the  cross 
flow  as  will  be  shown  later. 


III.  DRAG,  LATERAL  FORCES  AMD  MOMENTS 


In  order  to  determine  the  aerodynamic  forces  on  the  body, 
some  information  about  the  order  of  magnitude  of  some  pertinent 
quantities  is  needed.  Following  Ward's  proof,  we  remark  first 
that  on  the  body  the  normal  and  tangential  velocity  components 
in  a  cross  plane  (perpendicular  to  x-axis)  are  both  0(0, 


consequently  on  and  near  the  body 


0  CO 


then  Eqs.  (6)  and  (12)  indicate  that 


^  o  (  t*  )  )  n  =  0,  1,  2,  . . . 


ana  .  *2.  c 

TZ  '  ^  o^h[) 

a/  o  (  ^  ) 

^  (16) 

A  '^4  Art') 

'  a-  zb  "  0  C  } 

In  other  words,  the  order  of  magnitude  of  these  functions  are 


o  C  1  ) 


not  changed  due  to  the  nonequilibrium  dissociation. 


(a)  Pressure  relation 

In  deriving  Eq.  (1),  we  used  the  linearized  Euler's  equa¬ 


tion 


* 
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which  is  just  the  linear  approximation  to  the  pressure.  This 
pressure  relation  is  usually  good  for  two  dimensional  flow. 

Pnr  a  slendp/r  hndv  ««  up  have  ah  run  l-haf  anrf 


For  a  slender  body  as  we  have  just  shown  that  and  (— 2- 

1  Vo 

(5^)  are  of  the  same  order  of  magnitude  on  and  near  the 


body,  a  better  approximation  is  therefore 


-  \v 


a.  \4>  j_  /^<j> 


+  } (17) 


which  is  the  familiar  quadratic  approximation  of  the  pressure 
for  the  slender  body  ,  That  the  linearized  equation  (1),  where 
the  terms  like  )  etc.  are  neglected,  in  combination  with 

the  quadratic  approximation  of  Eq.  (17),  where  (^~)  and 

are  kept,  gives  the  correct  first  approximation  for 
the  aerodynamic  forces  on  the  body  is  well  known  in  the  conven¬ 
tional  slender  body  theory.  We  here  merely  demonstrated  that 
the  same  pressure  Eq.  (17)  should  be  used  for  the  same  reason, 
namely  the  terms  on  the  right  side  are  all  of  the  same  order  of 
magnitude  near  the  body. 


(b)  Density  relation 

In  deriving  Eq.  (1) ,  we  also  used  the  linearized  continuity 


equation 


u  if  LJs 


)  +  = 


meanwhile  Eqs.  (5)  and  (11)  indicate  that  in  our  present  non¬ 
equilibrium  case,  also  satisfies  the  Laplace's  equation  in 


% 


Che  cross  plane  i.e., 


hence 


4  1  jL£ 

'b  ^  3i  "ci  JJj, 


O 


*> 


=  o 


or 

-t  *,  ,  _  _L  >$ 

-v  'b^o  (18) 

Y * 

which  Is  again  Che  same  form  as  used  In  Che  convenclonal  case. 
Having  Chus  shown  Chac  Che  same  pressure  and  densicy  relacions 
as  in  Che  convenclonal  case  can  be  used  for  our  presenC  non¬ 
equilibrium  case,  we  are  now  in  a  position  Co  calculate  Che 
aerodynamic  forces. 


(c)  Drag  relation 

Our  solution  <£>  given  by  Eqs.  (5),  (11),  pressure  Eq.  (17) 
density  Eq.  (18)  are  of  the  identical  forms  as  these  for  the 
conventional  case,  it  follows  that  the  expressions  for  the  drag, 
lateral  forces  and  moments  deduced  by  Ward^  from  the  general 
momentum  theory  apply  equally  to  our  present  nonequilibrium  flow 
the  drag  is  given  by 


2 .T 

yr-'i. 


<*e(.0  &o( 0  +  A  t  fao'CO 

"^•6-  J 


(19) 
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where  C  denotes  the  contour  of  the  cross  section  of  the  body 
In  a  plane  normal  to  the  x>axis.  ^  denotes  the 

derivative  of  the  perturbation  potential  taken  with  respect 
to  the  outward  normal  to  the  body  cross  section  contour  (Fig.  2). 
is  the  base  pressure  coefficient. 


For  supersonic  flow,  substitution  }  4K>CO 

from  Eqs.  (7a),  (7b)  gives  under  the  assumption  of  pointed  nose 


i  ^  "'V 
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(20) 


where  the  first  four  terms  are  same  as  those  obtained  by  Ward. 
The  last  two  are  purely  due  to  the  nonequilibrium  dissociation. 
The  third  term  contains  the  contributions  due  to  the  angle  of 
attack  and  also  the  nonequilibrium  dissociation  because  <j> 

Qro  .  This  term  will  be  evaluated  later  when  a 


contains 


specific  type  of  body  like  the  body  of  revolution  is  to  be 
considered. 


18 


It  is  interesting  to  notice  that  for  pointed  nose  slender 
bodies ,  Eq.  (20)  does  not  depend  on  the  Mach  number,  a  fact 
which  is  well  known  in  the  conventional  case.  A  consequence 
of  this  Mach  number  independence  is  that  to  the  present  approxi¬ 
mation,  a  slender  body  experiences  the  same  drag  in  an  equili- 
brium  or  frozen  dissociating  gas  flow  as  in  the  conventional  gas 
flow,  because  the  governing  equations  for  these  cases  differ 
only  in  Mach  number. 


There  are  two  cases  for  which  Eq.  (20)  may  be  simplified, 
(a)  when  S(  [  )  =0  i.e.,  the  body  is  pointed  at  both  ends, 
and  (b)  when  S'(  £)  =0  and  the  generators  of  the  cylinder  are 
parallel  to  the  free  stream  which  makes  ^  -O.  For 

TV  CAS 

both  of  these  cases,  omitting  the  base  pressure  for  case  (b) , 


±  TT 
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(21) 


0  0 

Since  ^,  >  quantity  inside  the  bracket  in  the 

second  integral  is  always  positive.  However  the  sign  of  S'* (x) 
and  S’(x)  is,  in  general,  dependent  the  body  shape.  The  sign 


of  the  second  integral  therefore  cannot  be  asserted,  in  other 
words,  one  cannot  say  whether  it  represents  a  drag  or  a  thrust 
in  general. 


Though  one  can  integrate  the  second  integral  in  Eq.  (21) 


by  parts,  and  using  the  condition  S'(0)  8  0,  to  yield 

.  „  .r.  ..  .v  .  _>*  »?..*♦  .« 


L *  r  «  «3 

j  j  -  -j—f-  C*-*0  1 1L 


--  ]«M*(21a) 


or,  by  the  definition  of  exponential  integral,  one  can  rewrite 

Qg  |  J 

+|  J  A*)  .$"<£)  -A)}  -  d^cbc 


=^s'a)^h  4 1  ]]<**<* 


(21b) 


it  still  does  not  appear  to  be  possible  that  a  definite  con¬ 
clusion  on  the  sign  may  be  drawn-  from  any  one  of  these  alterna¬ 
tive  expressions  in  general.  The  situation  is  different  for 
the  first  integral  in  Eq.  (21)  due  to  the  different  behavior  of 
log  (x )  and^£l[-f“a(*---0]  -  This  first  integral 

certainly  gives  a  drag  and  is  larger  in  magnitude.  Some  further 
discussion  about  the  sign  based  on  the  mean  value  theorem  is 
given  in  the  Appendix. 


In  the  case  of  a  slender  cone,  R(x)  8  (constant)x. 

It  follows  then  that 

S'O)  *  2"iT  R  V  >  O 

SM  (V)  -  27T  V  ^  o 
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therefore  S11  Cx)  S*(4)  >  o  and  the  second  integral 

of  Eq.  (21)  represents  an  increase  in  drag. 


For  subsonic  flow,  we  consider  the  body  to  be  pointed 


at  both  ends  so  we  may  assume  the  absence  of  the  upstream 

influence  of  a  wake.  Eq.  (19)  thus  becomes 

_  i 

4  !<  r  . 
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Q-o  (*  )  ^0CO 


(19a) 


substituting  a  (x\  from  Eqs.  (7a),  (15), 


have 


rA»>s'(<) 


-^2^-4)  _  ^  (x  -J$  ) 


3C  -  4 


the  other  terms  either  vanish  by  the  assumption  of  pointed  ends 
or  cancel  out  each  other.  Eq.  (22)  says  that  there  exist  a  non¬ 
equilibrium  subsonic  drag  or  thrust  which  in  fact  has  the  same 
form  as  its  counterpart  in  supersonic  flow  given  in  Eq.  (21) 
except  that  the  sign  is  just  reversed,  because  for  subsonic  flow, 
(?*  4  B  a  the  quantity  inside  the  bracket  in  Eq.  (22)  is  always 
negative.  The  possibility  of  a  nonequilibrium  thrust  at  subsonic 
speeds  as  indicated  here  deserves  further  attention.  In  fact, 
as  pointed  out  in  Introduction,  Clarke^  has  obtained  a  subsonic 
nonequilibrium  drag  formula  same  as  in  Eq.  (22)  but  he  has  con¬ 
cluded  that  the  slender  body  would  experience  a  positive  non¬ 
equilibrium  drag.  To  interpret  the  present  result,  we  need  to 
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obtain  some  physical  insight  of  the  nonequilibrium  drag  problem. 
This  is  exactly  what  has  been  accomplished  in  Ref.  8.  We  shall 
discuss  the  physical  aspects  of  the  nonequilibrium  drag  probelm 
in  Section  IV. 

(d)  .  Lateral  forces  and  moments 

The  lateral  forces  Fz,  Fy  acting  on  the  body  can  also  be 
determined  from  the  momentum  considerations.  However  it  is 
more  convenient  to  write  them  in  complex  form 

F  ~  F-^  +  1  F2 

so  that  the  theory  of  complex  function  can  be  applied.  Clearly 
the  real  part  Fy  in  our  present  coordinate  system  represents 
the  lift.  The  expression  for  F  deduced  by  Ward  depends  only  on 
the  coefficient  (x )  and  the  base  cross  section,  but  not  on 

the  term  (jc)  .  Consequently  the  lateral  forces  and  their 
moments  about  the  y  and  z-axis  are  the  same  for  both  supersonic 
and  subsonic  conventional  gas  flow. 

For  our  present  nonequilibrium  flow,  the  nonequilibrium 
effects  as  discussed  above,  also  enter  only  through  the  term 
bQ(x),  Since  Ward's  expression  for  the  lateral  forces  and 
moments  apply  equally  to  our  case,  therefore  we  conclude  the 
nonequilibrium  dissociation  does  not  affect  the  lateral  forces 
and  moments.  To  emphasize  this  point  we  may  state  that  a 


slender  body  according  to  the  present  approximation  experiences 
the  same  lateral  forces  and  moments  regardless  of  the  flow  being 
supersonic  or  subsonic,  dissociating  or  non -dissociating,  in 
equilibrium  or  frozen  or  in  nonequilibrium.  This  conclusion  is 
not  surprising  because  the  cross  flow  for  all  these  cases  is 
regarded  as  incompressible  flow  in  the  present  approximation. 


23 


IV.  PHYSICAL  MECHANISM  OF  NONEQUILIBRIUM  DRAG 


Calculation  in  Section  III  show  that  a  slender  body  would 
experience  an  additional  drag  (positive  or  negative,  depending  on 
the  body  shape)  caused  by  chemical  relaxation  in  a  dissociating  gas. 
The  recent  study  of  Kusukawa  and  Li  provides  an  explanation  of  the 
physical  mechanism  of  this  nonequilibrium  drag.  They  showed  that 
the  body  drag  can  be  computed  as  follows: 


where 
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_i  - 
?  ~ 


5  Y  v 


F  5  b'1-' 
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(23) 


2  =  and  2  -  +  "W  denote  two  points  P  and  Gl  on  a  streamtube, 
upstream  and  downstream  far  from  the  body  where  all  perturbations 
almost  vanish  and  denotes  the  integration  with  respect  to 

all  streamtubes,  the  streamtube  cross  section  being  dUr  .  From 
Eq.  (23)  it  is  seen  that  the  sign  of  the  drag  is  determined  by  the 
sign  of  the  following  integral  along  a  streamtube: 

2=  -tV 


f  d  V  = 


2  =  +  v> 
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2--V 


( K.+ 1')  <Av  ' 


(23a) 
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where 


=  t  -  K 


v  -  V  -  vv 

The  subscript  y>  denotes  here  free  stream  quantities  and  the 
primed  quantities  are  perturbations*.  Therefore  if  we  can  determine 
^  and  y  /  along  a  streamtube  we  shall  be  able  to  compute  the  drag. 
Using  the  Lagrange's  approach,  Kusukawa  and  Li  obtained  a  linear 
equation  that  provides  the  governing  relation  between  the  variation 
of  ^  and  V’of  a  moving  material  element  on  a  streamtube.  In  the 
case  of  subsonic  flow  past  a  slender  body,  the  function  V'  between 
P  and  Ck  can  be  generally  represented  by  a  Fourier  series: 


V  - 


'T\ 


**  =  I 


s  v/v  ST- A  ~  +  £ 

*  Tl-> 


-tun 

of 


(24) 


where  t  is  the  time  variable,  ^  *  denotes  the  time  interval  during 
which  the  material  element  moves  from  P  to  &  .  The  Fourier  co¬ 

efficients  an  and  bn  must  be  determined  by  the  dynamic  and  kinematic 
conditions  appropriate  to  the  problem  and  indeed  they  must  be  de¬ 
pendent  on  the  body  shape.  In  Ref.  8,  it  has  been  shown  that  the 
perturbation  pressure  that  corresponds  to  a  Fourier  component 


V. 
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*  This  should  not  be  confused  with  the  prime  that  denotes 
differentiation  such  as  S'{  \)  in  Eq.  (7a). 


X 

From  Eqs.  (25)  and  (26),  we  conclude  that  \(/  and  are  out  of 


phase  due  to  relaxation  effects.  In  equilibrium  and  frozen  flows, 
however,  Vj/and  pTare  in  phase  and  =  C  .  The  consequence  of 
this  phase  shift  due  to  relaxation  effects  is  to  cause  a  nonequili¬ 
brium  drag  as  can  be  verified  from  application  of  these  results  in 
Eq.  (23a) .  Indeed,  in  equilibrium  and  frozen  flows,  =  °  >  ^  *  °  > 
and  in  nonequilibrium  flow,  j  O  ,  Kusukawa  and  Li  obtained 

ID  =  ^  Z  (  >°  (27) 

where 
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Therefore,  the  phase  shift  ^  for  the  functions  and  can  be 
regarded  the  physical  cause  of  the  appearance  of  a  nonequilibrium 
drag.  This  drag  is  positive  by  Eq.  (27)  because  an  and  bn,  which 


are  dependent  on  the  body  shape,  have  been  assumed  to  be  real  quan¬ 
tities.  This  shows  that  for  a  class  of  bodies,  for  which  y '  can 
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be  represented  as  in  Eq.  (24)  with  real  an  and  bn,  we  must  expect 
the  subsonic  nonequilibrium  drag  to  be  positive.  This  also  shows 

that  the  possibility  of  J)  3 r  exists  only  if 

1  ‘2- 

^  which  implies  that  an  and  bn  must  be  complex  quan¬ 

tities  for  such  bodies.  This  possibility  has  not  been  seriously 

/ 

considered  here  since  v  in  Eq.  (24)  is  a  real  physical  quantity 
and  generally  an,  bn  are  expected  to  be  real  quantities.  We  remark 
therefore  that  the  subsonic  nonequilibrium  drag  on  a  doubly  pointed 
slender  body  is  positive.  In  this  sense  the  integral  in  Eq.  (22) 
should  then  take  negative  sign.  Kusukawa  and  Li  did  not  consider  the 
supersonic  flow  case  in  detail.  They  pointed  out  that  in  such  a  case 
the  drag  would  consist  of  supersonic  wave  drag  and  the  non equilibrium 
drag.  The  nonequilibrium  effects  would  tend  to  provide  damping  of 
the  disturbances  propagating  along  Mach  waves  and  would  thus  be  ex¬ 
pected  to  decrease  the  wave  drag  from  the  conventional  value. 
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V.  APPLICATION  TO  BODIES  OF  REVOLUTION 

Since  the  lateral  forces  and  moments  are  same  as  the  conven¬ 
tional  case,  we  shall  consider  here  only  the  drag  problem.  Fur¬ 
thermore  for  subsonic  flow,  the  body  must  be  pointed  at  both  ends 
so  that  our  present  theory  may  be  applicable,  and  the  drag  has  been 
completely  determined  in  Eq.  (22).  In  *the  following  we  therefore 
consider  the  supersonic  flow  only,  and  the  body  may  have  a  flat  base. 

Our  method  of  solution  is  first  to  obtain  a  complex  potential 
W  -  Causing  our  knowledge  in  the  classical  incompressible  two 
dimensional  flow.  When  a  body  of  revolution  is  at  a  small  angle  of 
attack  ~«^0  (Fig.  Ill),  we  may  consider  the  cross  sections  in  a  plane 
normal  to  the  x-axis  to  be  circles  through  they  actually  are  ellipses 
with  an  error  of  O  .  It  is  well  known  that  for  a  circular 

cylinder,  the  complex  potential  should  be  the  superposition  of  a 
source  and  a  doublet,  i.e. 

w  -  £,(-0  =  (do  (28) 

*Tr  °  T  5-5, 

where  is  the  center  of  the  cross  section,  given  by  %0  -  _ 

denotes  the  angle  of  attack.  The  coefficients  were  so  chosen 

that  the  boundary  conditions  will  be  satisfied  as  shown  below. 
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On  the  body  Z>  -  substituting  SC*>  » 

Z>q-  -^0^c  in  Eq.  (28),  we  obtain 

W^.u.e)  =  t.eO  +  ^|^.)[^3TtCxHi.e]  +  ^fe)e  (29) 


The  real  part  in  Eq.  (29)  is 

<p  =  g0  oo  +  — -  dc  t^  t<  (-I  ^  +  i^Y> "RC '  ^  ^ 

From  the  boundary  condition,  we  have 
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for  a  body  of  revolution  at  angle  of  attack  °l0 


(30) 
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therefore,  we  obtain 

< p  _  ^  Q  f; 1  (  ...  _)  —  U  ,  c  C'\  9^j  (31) 

cm 

This  same  relation  can  be  obtained  by  differentiating  Eq.  (30). 
Thus,  the  coefficient  in  Eq.  (28)  have  been  chosen  correctly. 
Furthermore,  expanding  Eq.  (28) ,  we  obtain 
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which  shows  that 

a,  oO  =■  -  “U, 
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while  Cl  is  the  same  as  given  by  Eq.  (7a). 


Known  <f>  and  ^,,£-  ,  we  can  evaluate  the  contour  integral 


f  *p  |  *-n  t^le  drag  expression  Eq. 
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Substituting  (^  (  £  )  from  Eq.  (7b)  in  the  above  expression  gives 
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We  may  write  the  drag  coefficient  as  follows: 
lb 

C  ^  =  ~  ~  C-r>  +  Ct)  ■+■  C  r~. 

*  °  1  *'*«* 


(33) 


where  C>>  ,  Op  are  respectively  the  drag  coefficient 

■*'6  J  i  ■  "Y\  err  (3y 

at  zero  angle  of  attack,  the  induced  drag  coefficient  and  the  non¬ 


equilibrium  drag  coefficient. 
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Substitution  of  Eq.  (32)  into  Eq.  (20)  then  yields 

Cp.  =  -fT,  II  sV>  s"(0  h  \*-A\  ‘M31 

°  °  |  i 
i  s!il>  [  s'«)  ic a  < >-*, )  ^  - 52)  io3  ^  e«  ) 


Cp,  --  ^  s«) 


(34b) 


^™,e,  *  ^  (f5Vs'«) 

y  oo 


e-  cA^d 


X.-Z, 


_  s/l) 


s&) 

2U  j 

o 


e-&04) 


(34c) 


-  e 


A  ^ 


We  have  thus  demonstrated  that  the  drag  contributions  due  to  the 
angle  of  attack  and  the  nonequilibrium  dissociation  can  be  separated 
as  might  have  been  expected  in  the  present  linearized  theory. 
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VI  Conclusions 

In  the  present  paper  we  have  considered  the  inviscid  com¬ 
pressible  linearized  dissociating  gas  flow  over  slender  bodies 
of  general  cross  section.  Essentially  it  is  an  extension  of  the 
earlier  works  of  Ward  and  Sears  and  Adams  to  include  the  disso¬ 
ciation  effects.  The  transform  methods  have  been  used  to  obtain 
the  solutions,  both  supersonic  and  subsonic  case  have  been  examined. 

The  effects  of  Mach  number  and  the  rate  of  dissociation  enter 
only  through  the  term  $-0(;x).  The  pressure  and  density  relations 
have  been  shown  to  take  the  same  forms  as  in  the  conventional 
theory  (c.f.  Eqs.  (17),  (18)  ),  consequently,  the  expressions  for 
the  lateral  forces  mements  and  drag  deduced  by  Ward  from  general 
momentum  theory  apply  equally  to  the  present  case. 

For  the  equilibrium  and  frozen  cases,  the  lateral  forces  moments 
and  drag  were  found  identical  to  those  for  the  conventional  case  pro¬ 
vided  that  the  tail  end  of  the  body  is  either  pointed  or  has  a  cylin¬ 
drical  form. 

For  the  nonequilibrium  case,  the  lateral  forces  and  moments 
were  found  to  be  the  same  as  in  the  conventional  case,  but  there  exist 
additional  nonequilibrium  terms  in  the  drag  expression.  They  may 
represent  drag  or  thrust.  If  the  body  is  pointed  at  both  ends,  the 
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nonequilibrium  term  is  the  sam  for  both  supersonic  and  subsonic 
flows  except  that  the  sign  is  just  opposite.  By  considering  the 
physical  mechanism  of  nonequilibrium  drag,  we  obtain  that  in  the 
subsonic  case,  it  represents  a  drag,  in  the  supersonic  case  it  is 
a  thrust. 

Application  to  bodies  of  revolution  in  the  case  of  supersonic 
flow  has  been  carried  out  in  detail.  The  drag  expression  can  be 
separated  into  three  parts,  (1)  the  drag  at  zero  angle  of  attack, 
(2)  the  induced  drag  and  (3)  the  nonequilibrium  terms  which  may 
represent  drag  or  thrust. 
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Appendix 


Estimation-  of  Integral  by  Mean  Value  Theorem 


Mean  Value  Theorem 

Let  (a)  and  (x)  be  bounded  and  integrable  in  £a,bj 

and  (b)  ^  (jx)  keep  the  same  sign  throughout  j_a,bj,  then 

|  f  (*)  cLx  ^  j  ^  (A  J  q(x)dx 
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where  M  and  m  are  the  upper  and  lower  bounds  of  f('Oin  [a,bj  , 


or  we  write 


S  l 

|  f(x)cl'x  =  y(y)  <^(x)dx  ) 


For  proof  of  this  theorem,  see,  for  example,  Courant's  Differential 
and  Integral  Calculus,  Vol.  I,  p.  127,  2nd  ed. 

We  now  apply  this  theorem  for  the  estimation  of  the  integral 
inEq.  (22).  We  know  for  the  subsonic  flow  ^ (  <  0^ 
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Now  since 


[  E  »  C“ ^2  i  l  (-  3.)  +  j  >  O 


we  may  again  apply  the  mean  value  theorem  so  that 

^Iccj^Ddx  0^*{ 
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The  last  integral  is  positive,  but  both  5^6)  and  S/ (j<0')may  be 

positive  and  negative  depending  on  the  body's  shape,  consequently 
T)  may  be  positive  or  negative,  i.e.  drag  or  thrust. 
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